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Abstract 

The problem of reconstructing functions from their asymptotic expansions in pow- 
ers of a small variable is addressed by deriving a novel type of approximants. The 
derivation is based on the self-similar approximation theory, which presents the pas- 
sage from one approximant to another as the motion realized by a dynamical system 
with the property of group self-similarity. The derived approximants, because of their 
form, are named the self-similar factor approximants. These complement the obtained 
earlier self-similar exponential approximants and self-similar root approximants. The 
specific feature of the self-similar factor approximants is that their control functions, 
providing convergence of the computational algorithm, are completely defined from 
the accuracy-through-order conditions. These approximants contain the Pade approx- 
imants as a particular case, and in some limit they can be reduced to the self-similar 
exponential approximants previously introduced by two of us. It is proved that the 
self-similar factor approximants are able to reproduce exactly a wide class of functions 
which include a variety of transcendental functions. For other functions, not pertain- 
ing to this exactly reproducible class, the factor approximants provide very accurate 
approximations, whose accuracy surpasses significantly that of the most accurate Pade 
approximants. This is illustrated by a number of examples showing the generality and 
accuracy of the factor approximants even when conventional techniques meet serious 
difficulties. 



1 Introduction 



The problem of reconstructing functions from their perturbative asymptotic expansions in 
powers of a parameter or a variable is so frequently met in physics and in applied sciences 
that there is no necessity to explain its importance. The best known methods for such a 
reconstruction are the Pade approximation and Borel summation, including their variants 
and combinations [1,2]. These techniques usually require that a large number of terms of 
an asymptotic expansion be available. The Borel summation demands, in addition, that 
the high-order expansion coefficients be given and the analytic properties of the sought 
function on the complex plane be prescribed. However, the overwhelming majority of realistic 
physical problems are too complicated and perturbation theory is only able to derive a few 
first terms. And the luxury of knowing in advance the analytic properties of an unknown 
function, together with its high-order expansion coefficients, as is required for the Borel 
summation, is practically never available. Because of the latter, Pade approximants are more 
often employed in applications, although their usage also confronts with several difficulties, 
among which the most notorious are the appearance of spurious poles and the poor recovery 
of non-integer critical exponents. 

An alternative approach to the problem of reconstructing functions has been developed, 
whose basic ideas are as follows. First of all, to improve the convergence property of a per- 
turbative sequence, it is necessary to introduce control functions defined by an optimization 
procedure [3-5] . This idea makes the foundation of the optimized perturbation theory that is 
now widely employed for various applications [3-16]. The second pivotal idea is to consider 
the successive passage from one approximation to the next one as a dynamical evolution 
on the manifold of approximants, which is formalized by the notion of group self- similarity 
[17-22]. And the third principal point is the introduction of control functions in the course 
of rearranging perturbative asymptotic expansions by means of algebraic transforms [23-27] . 
Because of their specific scaling properties typical of fractals [28,29], the algebraic transforms 
can also be called fractal transforms [30]. By using this technique, two types of approximants 
have been obtained, self-similar exponentials and self-similar roots [24-27] . 

In the present paper, we suggest a different approximation scheme resulting in what may 
be named self-similar factor approximants. These new approximants possess an important 
principal property distinguishing them from the self-similar approximants mentioned above: 
the control parameters, entering the self-similar factors, can be completely defined from a 
given asymptotic expansion by the so-called "accuracy-through-order" matching method. 
This is in contrast with the self-similar exponentials whose controls, designed to improve 
convergence, are defined from additional optimization conditions. This method of accuracy- 
through-order is also different from the determination of the control parameters of the self- 
similar roots which are determined by matching two asymptotic expansions valid in the 
neighborhood of two different asymptotic points. Being based on the sole initial asymptotic 
expansion, the self-similar factors have the advantage of simplicity which makes their usage 
quite attractive. Furthermore, these approximants allow one to reconstruct exactly a wide 
class of functions. And when they do not yield exact answers, they provide very accurate 
approximations, essentially more accurate than given by Pade approximants. 

We, first, give in the next section the mathematical foundation of the self-similar factor 
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approximants. Sections 3 and 4 are devoted to examples chosen both for their illustrative 
properties, the difficulties they pose to the more convential Pade technique and the relevance 
to several physical problems (polymers, state- and velocity-dependent solid friction dynamics, 
critical phenomena in field theory and in Ising models). 

2 Mathematical Foundation 
2.1 Derivation of factor approximants 

Assume that we are solving a complicated problem, aiming at finding a real function f{x) 
of a real variable x. Because of the complexity of the problem, the only thing we are able 
to do is to invoke a kind of perturbation theory for obtaining approximate expressions fk{x) 
of order k = 0,1,2,..., valid in the asymptotic vicinity of a; = 0. Usually, the asymptotic 
approximants fk{x) can be presented as a power series of x and written in the form 

k 

fk{x) = /o ^k{x) , (Pk{x) = X] ' (1) 

n=0 

where ipk{x) is a dimensionless function with Oq = 1. Writing down the sought function as 
an asymptotic series 

/(x)~/oX:a„x" + ... , (2) 

n=0 

where x ^ 0, gives little consolation, since in real problems x is rarely asymptotically small, 
but usually it is finite and may be even very large. How could we reconstruct the function 

f{x) for finite values of x from the only knowledge of its asymptotic expansion? 

An answer to this question can be provided by the self-similar approximation theory 
[17-22], with control parameters introduced by means of the fractal transform [23-27,30], 
defined as 

Fk{x,s) = x' fk{x) . (3) 

This leads to the self-similar exponential and self-similar root approximants. Now, we shall 
follow a slightly different procedure, which is actually motivated by the very idea of group 
self-similarity underlying the construction of self-similar approximants. 

For a more efficient use of the group self-similarity, we propose to present an initially 
given asymptotic expansion in the most symmetric way. To this end, we introduce the 
factor functions 

(Pkp{x) = 1 + bkp X (p= 1,2, A; > 1) . (4) 

Let us consider the finite series (1) as a polynomial over the field of complex numbers. Then, 
by the fundamental theorem of algebra [31], such a polynomial can be split in a unique way 
into a product of the irreducible factors (4) , so that 

k 

Vk{x) = n ^kp{x) (/c = 1, 2, . . .) . (5) 

p=l 
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The representation (5) possesses the scahng property, for which if ip^p ^^kp, then ipk 
X^fk- Such a scahng property is the simplest rudimental example of functional self-similarity. 
In this way, the sum (2) can be rewritten as a product 

k 

fix) ~ /o n i'Pkpix) + ...]. (6) 

p=l 

Now, instead of accomplishing the self-similar summation for the whole right-hand side of 
Eq. (2), we may perform it for each factor in the product (6). Thus, we define the fractal 
transform 

^kp{x, s) = x" (pkp{x) , (7) 

construct an approximate cascade whose trajectory is bijcctivc to the sequence {^^^(x, s)}, 
embed the cascade into an approximation flow, integrate the flow evolution equation, and 
realize the inverse fractal transform. All this machinery, with all details, has been expounded 
in previous papers [17-27], and we therefore do not repeat it here. As a result, each factor 
(4) can be shown to be renormalized into 

^l^{x) = (1 + Akp xr^ , (8) 

where A^p and Ukp are control parameters, or simply controls. And a /c-order approximation 
for the sought function f{x) is given by the self- similar factor approximant 

rki^) = foU^Ux). (9) 

p=i 

The controls A^p and Ukp are determined by expanding the approximant (9) in powers of x 
and comparing this expansion with the series (2). For short, this can be called a re-expansion 
procedure, which sometimes is also named the accuracy-through-order relationship. The 
equations defining the amplitudes A^p and exponents rikp can be cast in the form 

j^rikpAl^i-ir^'nb^, (10) 
p=i 

where 

6„=- hm — In^a^x™ . (11) 

As is easy to check, Eqs. (10) and (11) follow from equating the asymptotic expansions for 
the logarithms of the factor approximant (9) and of series (2). For each given k, there are 
2k unknowns in the left-hand side of Eq. (10). Hence, in the series (2), one should have 
2k nontrivial terms, which makes n = 1,2, ... ,2k. The series of odd orders 2k + 1 can 
also be processed, for which one needs to consider f{x) — /o, instead of f{x). The factor 
approximants, based on even and odd numbers of terms of a series with alternating signs, 
often bracket the sought function from below and above. This bracketing is analogous to 
that occurring for self-similar exponential approximants based on even or odd numbers of 
asymptotic terms [26,32]. 
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The control parameters may be complex-valued, since, for obtaining the factorized form 
(5), the sum (1) was treated as a polynomial over the field of complex numbers. But, since 
the considered function is real, all complex-valued factors should arise in complex conjugate 
pairs, so that their product be always real. 

2.2 Exactly reproducible class of functions 

The structure of the self-similar factor approximants (9) suggests that there exists a whole 
class of functions that are exactly reproducible by means of these approximants. This class 
is defined as follows. Let P„(x) be an irreducible polynomial in a real variable x of degree 
n over the field of real numbers and let and (3j be complex numbers. Compose the 
real- valued products of powers of such irreducible polynomials as Jlj ^Mii^) ^-^id Qnj{^)j 
where J2i = ^ ^ind J2j = This implies that complex powers, if any, always come 
in complex conjugate pairs. Let these products have no common divisors, such that the ratio 

be irreducible. Denote by 7^ a class of functions, which is composed of all products of the 
forms (12) that play the role of the prime representatives for this class. 

Theorem. A function f{x) can be exactly reproduced by the self-similar factor approxi- 
mants (9) if and only if this function belongs to the class TZ, with the prime representatives 
(12) being exactly reproducible by /^(x) in all orders k > M + N. 

Proof. To prove the proposition for the whole class TZ, it is necessary and sufficient to 
prove it for the prime representatives (12). According to the fundamental theorem of algebra 
[31], each polynomial in one real variable, over the field of real numbers, can be split into 
factors of the first and second degree, over the field of real numbers, and into factors of 
the first degree, over the field of complex numbers. This allows us to split each polynomial 
entering the prime representative (12) into the product of first-degree factors with complex 
coefficients. Since, by definition, the ratio (12) is irreducible, it can be presented as a 
product of M + first-degree factors. After this, it acquires the form identical to the factor 
approximant (9) of the order k = M + N, taking into account that the exponents ra^p in 
(H) can be negative. Hence the latter pertains to the class TZ. And if a function is exactly 
reproducible by an approximant (9), this function must be from the class 71. 

By construction, all parameters (coefficients and powers) of the approximant (9) are 
defined by equating its asymptotic, in x — 0, expansion with that of the function (12). 
An asymptotic expansion in the sense of Poincare is uniquely defined by the function itself. 
Two functions, having an identical dependence on the variable and coinciding asymptotic 
expansions, coincide. That is, a function from the class TZ is exactly reproducible by a factor 
approximant (9). Finally, if a function f{x) exactly equals an approximant f^{x) of the 
order k = p, then the higher approximants, with k > p, derived from the same asymptotic 
expansion of the same function f{x), will coincide with each other. 

Remarks. The class TZ of the functions, exactly reproducible by means of the self-similar 
factor approximants, is significantly wider than the class of rational functions that can be 
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exactly reproduced by Pade approximants. In addition to rational functions, the class TZ 
also includes transcendental functions. Because of this, the self-similar factor approximants 
should provide a better accuracy for a wider class of functions, as compared to Pade approx- 
imants. In what follows, we shall illustrate this by a variety of examples. 

For a given real- valued asymptotic series (2) of a real function f{x), the factor approxi- 
mants (9), by construction, are real in the asymptotic region of x — > 0. However, they may 
become complex for finite x. If an approximant f^ix) becomes complex for x > Xk, this 
implies that the region of validity of fk{x) is restricted by the interval [0, Xk\- 

2.3 Relation to the Park method 

Factors (8) are appropriate for describing the behavior of functions in the vicinity of critical 
points. A value x = Xc is termed a critical point of a function f{x) if at this point the latter 
is either zero, f{xc) = 0, or possesses an algebraic singularity, that is f{x) ~ (1 — x/xc)~^ as 
X — > Xci where /3 is positive. This fact was, actually, employed by Park [33] who suggested 
a method for locating the critical points and calculating the critical exponents. His method 
is formulated as follows. Assume that: (i) a real function f{x) has a critical point Xc, (ii) in 
the neighbourhood of the critical point, the function can be represented as 

f{x)^Ujl{l-Bpx)-^- {x^x,) , 
p 

with all Bp and I3p being real; (iii) the physical critical point corresponds to that which is 
the closest to the origin, such that, arranging Bp in the descending order of their absolute 
values, \Bp\ > one has Xc = -Bf^. Then, defining the coefficients 6„ by the expansion 

oo 

In fix) = J2bn x^, 

n=l 

one obtains 

h h"- 
Bi = hm , /3i = lim n . 

n— ►oo h , n— »oo rP~^ 

The proof of this statement is based on the generalized Polya theorem [34], which extends the 
theorem by Polya [35], initially proved for entire functions of genus zero, when all /3p = — 1, 
to the case of real j3p. The Park method for defining the critical points and critical indices is 
closely related to the Pade analysis of logarithmic derivatives of a series, though in the latter 
case there is no prior knowledge for the convergence of estimates from Pade approximants 
[36] , which would be analogous to the generalized Polya theorem. 

The principal difference of our approach from the Park method is in the following. First 
of all, we never require that the sought function be exactly factorizable, but we derive the 
form (9) as an approximation to this function. Second, we do not impose a constraint that 
the function must necessarily possess a critical point, and if so, we consider the function 
not solely in the neighbourhood of the latter, but in the whole region [0, Xc]. Third, since 
we deal with a much more general case, the amplitudes A^p and exponentials n^p are not 
compulsory real, but may be complex valued. 
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3 Examples 



Any function from the class TZ can be reproduced, according to Theorem, exactly, provided 
that there are enough terms in series (2). However, a reproducible function may be not 
exactly reproduced when the asymptotic expansion (2) contains not enough terms. In other 
words, a function to be recovered may be exactly reproducible in principle, but in practice, 
we may have access to only a few terms in the asymptotic expansion. Then an important 
question to ask is how well the factor approximants are able to approximate such a function 
and also it is interesting to observe how the factor approximants converge to the exact result. 
This problem will be considered in the next subsection. 

Those functions that are not from the class TZ cannot be reproduced exactly, but they 
can be very well approximated by the self-similar factor approximants (9), as we illustrate 
in the following subsections. 

3.1 Convergence to exact result 

Consider the function 

fix) = (1 + 2x)^/''il + + 0.5x)^^\l + O.lxy/^ , (13) 

which is from the class TZ. Its expansion of eighth order has the coefficients /o = 1 and 
oi = 3.692 , 02 = 3.521 , og = 0.410 , 04 = 0.025 , 

as = -0.091 , ae = 0.145 , 0,7 = -0.220 , as = 0.335 , 

whose behaviour is rather irregular. If we take into account only four terms of series (2), 
then we get the approximant /I (a;), with 

A21 = 1.986 A22 = 0.721 , 7121 = 1.562 , 7122 = 0.818 . 

And for the factor approximant f^{x), constructed by means of series (2) of sixth order, we 
find 

^31 = 2 , = 0.960 , A33 = 0.321 , 

7131 = 1.503 , 7l32 = 0.583 , 7l33 = 0.398 . 

The best Pade approximant that can be built of the sixth-order series (2) is -P[4/2] (x) , whose 
accuracy is compared with that of f^ix) in Fig. 1. As is evident, the factor approximant is 
essentially more accurate, although is not yet exact. But in the next order, we obtain fl{x), 
with 

^41 = 2 , ^42 = 1 , ^43 = 0.5 , ^44 = 0.1 , 

3 11 1 

Uai = -, ^42 = - , ^43 = - , 7144 = ^ , 

which coincides with the exact function (13). 
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3.2 Combination of functions from TZ with exponentials 

The combination of functions from TZ and exponentials are approximated with a very good 
accuracy. As an example, let us consider 

The choice of the coefficients and the power m is not important, since the factor multiplying 
the exponential pertains to the class TZ of exactly reproducible functions. For concreteness, 
let us take A — 0.5, B — 0, C — 0.5, D — 0.1, and m — 0.5. Then, in the asymptotic series 
(2), we have fo — 1 and 

oi = -1.750 , 02 = 2.419 , as = -3.659 , 

a4 = 6.060 , as = -10.499 , ag = 18.622 . 

The first-order approximant (9) is, clearly, too simple for providing a good approximation 
for complicated functions. Therefore, here and in what follows, we start the analysis with 
the second-order approximant. For the considered case, we get in the second order 

A21 = 1.976 , ^22 = -0.077 , 7121 = -0.471 , 7122 = 10.651 

and in the third order 

^31 = 1-945 , ^32 = 0.501 , A33 = 0.000986 , 

ngi = -0.500 , 7132 = 0.497 , 7133 = -1039 . 

Both approximants f2{x) and f^{x) perfectly reproduce function (14). 

For comparison, we construct the Fade approximants based on the same number of 
terms in the series (2). Among all possible Fade approximants P[m/n]{x)-i we select the most 
accurate for the case considered. Note that the best Fade approximants are not necessarily 
diagonal. Here, these are P[i/s\ {x) and -P[2/4] {x). The percentage errors of these approximants, 
together with the error of f^{x) are shown in Fig. 2. One can observe that the accuracy of 
/g (x) is incomparably higher than that of the best Fade approximants, whose errors grow 
fast with X, reaching amplitudes of the order of 100%. In addition, the approximant P[2/i]{x) 
becomes negative for x > 5, which is qualitatively wrong for the positive function (14). 

3.3 Exponential multiplied by functions not from the class 71 

As an example of a function having no factors from the class TZ, let us consider 

f{x) = tanh(a;) exp(— x) . (15) 
In its asymptotic series (2), we have /o = 1 and 

- 1 - 1 - 1 

6 6 
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For the second-order approximant (9), we find 

A21 = -0.350 + 0.587 i , ^22 = ^21 > ^21 = 0.036 + 0.831 i , 7122 = n*i ■ 

And for tlie tliird order, we obtain 

A31 = -0.00337 + 0.650 i, ^32 = A^i , ^33 = 0.071 , 

ngi = -0.871 + 0.023 i , 7132 = n^i , 7133 = -13.777 . 

Again, wc compare f^{x) witli the most accurate Fade approximants that for this case are 
P[i/5](a;) and P[2/4]{x). The corresponding percentage errors are shown in Fig. 3. Again, 
we see that has much higher accuracy than the best Fade approximants. The most 

accurate of the latter, P[i/5](a;), becomes negative for x > 6, which is quahtatively wrong. 
A combination of a logarithm and an exponential yields a non-monotone function 

f{x) = ln(l + x) exp(-x) . (16) 

The coefficients of the related asymptotic expansion (2) are /o = 1 and 

3 4 
oi = - 2 ' = 3 > 03 = -1 , 

_ 89 _ 83 _ 593 

- 120 ' - " 144 ' - 1260 • 

For the parameters of the factor approximant we get 

A21 = 0.930 , ^22 = 0.013 , 7121 = -0.466 , 7122 = -83.334 , 

and for those of /| (x) , we find 

A31 = 0.973 , A32 = 0.574 , A33 = 0.00212 , 

n3i = -0.363 , 7132 = -0.215 , 7133 = -483.556 . 

The most accurate Fade approximant here is P[i/5](a;). The corresponding percentage errors 
arc presented in Fig. 4, from where it is seen that f^{x) is significantly more accurate than 
^[1/5] (a;)- 

3.4 Functions not from the class 71 which converge to a constant 
at infinity 

The previous functions converge to zero at infinity. Let us now consider the function 
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which increases at infinity to a finite value. In the expansion (2), we have fo — 1 and 

ai = 0.5 , 02 = -0.25 , 03 = 0.146 , 04 = -0.091 , 

05 = 0.059 , oe = -0.038 , 07 = 0.025 , Og = -0.016 . 
For the second-order factor approximant, we find 

A21 = 0.570 , A22 = 1.097 , 7121 = -0.671 , 7122 = 0.805 , 

and for f^{x), we have 

^31 = 1.041 , A32 = 0.794 , A33 = 0.265 , 

7131 = 1.243 , 7232 = -0.922 , 7l33 = -0.235 . 

The best Pade approximant P[2/2]{x) is less accurate than f^ix), as is shown in Fig. 5. The 
accuracy becomes even better for /I (a;), for which 

^41 = 0.151 , ^42 = 1-023 , ^43 = 0.881 , Au = 0.516 , 

7141 = -0.150 , 7142 = 1.675 , 7143 = -1.195 7144 = -0.267 . 

Notice that, passing to higher approximations, the sum Y^pf^kp decreases, 

n21 + ^22 = 0.133 , 7l3i + 7132 + 7133 = 0.086 , 

7l4i + 7l42 + 7l43 + 7l44 = 0.063 , 

which is the correct trend, since in the limit it should be 

k 

hm i^kp^Q ■ 
''-^^ p=i 

The accuracy of the factor approximants can also be improved by assuming the validity of 
the condition J2p ^kp = for finite k. 

3.5 Classical example of 0^ theory with strongly divergent series 

A classical example of a strongly divergent series is provided by the asymptotic expansion, 
in powers of a coupling parameter, of generating functionals in field theory or of partition 
functions in statistical mechanics. The generic structure of such divergent expansions is 
exemplified by expanding, in powers of the coupling the generating functional 

I{g) = ^ / exp (-(^2 _ ^^4\ ^ (^g) 

typical of ifr" field theory. In the asymptotic series (2), with g instead of x, the coefficients 
are /o = 1 and 

(-1)" ^ / r 
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The latter grows with n — cxd as ~ n". Such a series is divergent for any finite g. For 
the second-order approximant ^{g), we find 

A21 = 19.141 , A22 = 4.859 , 7121 = -0.00862 , 7122 = -0.120 . 

And for mg), we obtain 

^1 = 31.220 , A32 = 13.317 , A33 = 3.464 , 

ngi = -0.000526 , 7x32 = -0.022 , 7133 = -0.125 . 

The most accurate Pade approximant, constructed with the same number of the asymptotic 
terms, is P[3/3](5'). This is compared with I^{g) in Fig. 6. 

Again, we see that the accuracy of the factor approximant highly surpasses that of 
the Fade approximant. This is especially noticeable for large couphngs g, where the Fade 
approximant completely fails. Really, P[3/3](5f) is finite for g —>■ 00, which is in contradiction 
with the behavior oi I{g) at large g, where the integral (17) tends to zero, 

I{g)- 1.023 g-'-^' (^^00). (19) 

The factor approximants also decrease in the strong-coupling limit as 

r2{g)c^ 0.806 g-'-''\ i;{g) 0.807 g'^-''' (^^00). (20) 

Let us emphasize that all parameters of the self-similar factor approximants are defined 
only through the weak-coupling expansion. And it looks almost mysterious that they can 
reasonably extrapolate the behavior at the strong-coupling limit. The accuracy of the ap- 
proximants (9) increases with their order. Thus, for the approximant I^ig), we arrive at the 
parameters 

A41 = 43.965 , A42 = 23.064 , ^43 = 10.294 , ^44 = 2.677 , 

7141 = -0.000025 , 7142 = -0.00259 , 7143 = -0.035 , 7144 = -0.124 . 
In the strong-coupling limit, this yields 

Il{g)c^ 0.810 g-'-''' (5^00). (21) 

Note that the diagonal Fade approximants are always finite at infinity, thus allowing the 
approximation of only a very narrow class of functions [37]. Contrary to this, the factor 
approximants immediately catch the correct behavior at infinity; their accuracy increasing 
with the approximation order. Hence, the self-similar factor approximants, being based on 
an asymptotic expansion at zero, can correctly reproduce the behavior of the sought function 
for the whole range of its variable, including the behavior at infinity. 
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3.6 Expansion factor of a three-dimensional polymer chain 

As another physical example, let us consider the expansion factor a{z) for a three-dimensional 
polymer chain with excluded-volume interaction, where z is a dimensionless coupling param- 
eter [38,39]. An asymptotic series of the form (2), derived by means of perturbation theory 
[38], yields the coefficients /o = 1 and 

ai = ^ , as = -2.075385396 , 03 = 6.296879676 , 

04 = -25.05725072 , 05 = 116.134785 , ag = -594.71663 . 

In the strong-coupling limit, it has been established numerically [37,38] that a{z) can be 
accurately represented by 

a{z) ~ 1.531 + 0.184 ^-^'^^^^ {z ^ 00) . (22) 

For the approximant «2(-2)) we have 

A21 = 6.064 , A22 = 2.962 , riai = 0.105 , 7122 = 0.235 , 

with the strong-coupling limit 

~ 1.560;2°-3^° + 0.151^-'^-^6° (^^00) (23) 

being very close to the exact numerical value. In the next order, the parameters of a^iz) are 

^31 = 7.019 , - 4.635 , ^133 = 2.262 , 

nsi = 0.033 , 7132 = 0.164 , 7133 = 0.151 . 
And the strong-coupling limit gives 

aliz) ~ 1.551 + 0.166 ^-^'^^^ (^ ^ 00) . (24) 

Both the coefficients as well as the powers of the strong-coupling divergence are very accurate, 
as compared to the numerically found behavior of a{z), with the percentage error of 1%. 
The best Pade approximant here is P[3/3](2;). In Fig. 7, the percentage errors of al{z) and 
P[3/3]{^) shown, compared with the numerically fitted [39] equation 

0.1172 



a(^) = (^1 + 7.524 ;2 + 11.06 ■ . (25) 



As is evident, the factor approximant dramatically outperforms the Pade approximant. The 
latter completely fails at large z 00, where it is finite, while a{z) diverges. To the 
contrary, the factor approximant a^{z) possesses the correct behavior in the strong-coupling 
limit. Moreover, not only the main asymptotic, as 2; — > 00, term is very accurate, but the 
next term, describing the so-called correction to scahng, is also of good accuracy. Of special 
interest is the strong-coupling exponent 

1 1 \na(z) . ^, 
!/=- + - lim , 26 

2 4 ^-00 \nz ^ ^ 
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which is a kind of a critical index for the polymer chain [40]. According to the Muthukumar 
and Nickel [38,39] numerical estimate, vmn — 0.5886. For the self-similar factor approxi- 
mants, we have 

= 0.585 , z/* = 0.587 . 

Recent numerical estimates [40,41] give the value v = 0.5877 ± 0.0006, which is very close 
to i/3. 



3.7 Nonlinear differential equation for state-dependent solid fric- 
tion 

Let us now demonstrate how the self-similar factor approximants can be employed for solv- 
ing nonlinear differential equations. Let us consider the Ruina-Dieterich law of solid friction 
between two sohd surfaces shding against each other (see [42,43]). The Ruina-Dieterich 
law involves the so-called state variable denoted here as / in a dimensionlcss form, which 
is usually thought to quantify the true area of contacts of the asperities of two solid sur- 
faces. The state variable / obeys the following simple non-linear differential equation (put 
in dimensionless form both for / and t) 

f = /3 - f-"" , (27) 

where f3 and m are parameters. Equation (27) possesses two qualitatively different types of 
solutions corresponding to m < 1 and m > 1. 

Consider, first, the case m < 1. Let m = 0.85, /? = 0.526, and the initial condition 
/(O) = 0.5. Equation (27) allows us to derive the short-time expansion for f{t), presented 
as the series (2) in powers of time t, simply by using the Taylor expansion formula and by 
taking successive derivatives of (27). This gives /o = /(O) and 

oi = -0.750 , as = 0.102 , as = 0.012 , 

04 = 0.00226 , a5 = 0.000455 , = 0.0000861 . 

By the standard procedure, we obtain the factor approximant /|(i), with 

A21 = -0.391 + 0.116 i , A22 = A*2^ , nai = 0.682 + 0.932 i , 7122 = n^^ . 

And the parameters for /|(i) are 

>l3i = -0.430 + 0.118 i, ^32 = A;^ , /I33 = -0.149 , 

7131 = 0.790 + 0.471 i , 7132 = 71^1 , 7l33 = -0.272 . 

The accuracy of the approximants can be checked against the direct numerical solution of 
Eq. (27). In Fig. 8, we compare the percentage errors of f^{t) with those of the best Fade 
approximants P[2/2]{t) and P[3/3](t). As we see, the factor approximant is considerably more 
accurate. 
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4 Critical Behavior 



Solutions to physical problems often display a critical behavior, when a function f{x) tends, 
at a critical point Xc, either to zero or to infinity. Self-similar factor approximants make it 
possible to describe these two types of critical behavior, providing accurate estimates for the 
critical points as well as for critical indices. If one is interested solely in the critical behavior, 
then our approach reduces to the Park method, as is discussed in Sec. 11. However, we would 
like to stress that the factor approximants not solely describe well the critical neighbourhood, 
but provide an accurate approximation in the whole region [0,a;c]. It is also important to 
realize that, in order to achieve an accurate description, one does not necessarily require the 
knowledge of a large number of terms in an asymptotic series, but only a few terms are often 
sufficient. This applies to both types of critical points, either zeros or singularities. 

4.1 Critical behavior near zero 

An example of the first type is the function 

fix) = {cos ^/^y^' , (28) 

tending to zero at the critical point Xc = 7r^/4 = 2.467 with the critical index 1/3. the 
coefficients of the asymptotic series (2) are /o = 1 and 

oi = -0.167 , 02 = -0.014 , 03 = -0.00355 , 

04 = -0.000982 , 05 = -0.000295 , oe = -0.0000937 . 

Notice that all coefficients here are negative. Such asymptotic series with all coefficients of 
the same sign arc known to be very difficult for any kind of resummation procedure [2]. But, 
with the factor approximants (9), there is no problem in approximating the sought function 
in the whole region from x — to the critical point. For the approximant we get 

A21 = -0.405 , A22 = -0.024 , 7121 = 0.334 , 7122 = 1.333 . 

In this approximation, the critical point Xc — \A2i\~^ — 2.469 is already close to the exact 
value 7r^/4 = 2.467. The same quality of results is obtained for the critical index ^21, which 
almost coincides with the exact index 1/3. For the next approximation we find 

A31 = -0.405 , ^32 = -0.044 , A33 = -0.00541 , 

nsi = 0.333 , n32 = 0.375 , 7133 = 2.791 . 

The approximant f^{x) is practically indistinguishable from function (28) in the whole region 
[0,Xc]. The critical point Xc = |^3i|~^ = 2.467 yields the exact value 7r^/4 with a precision 
of about 10~^. The corresponding critical index 7131 differs from the exact index in the 
seventh decimal digit. The best Fade approximant P[3/2](x) is much less accurate, giving the 
critical point Xc = 2.507 and a very inaccurate critical index equal to 1 instead of 1/3. The 
accuracies of f^{x) and of P[3/2](a;) are compared in Fig. 9. 
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4.2 Critical singular behavior 

In several physical problems [42], the critical behavior is described by the function 

which diverges at the critical point Xc = 1, with the critical index 1/2. The asymptotic series 
(2) has the coefficients /o = 1 and 

ai = 0.637 , a2 = 0.405 , ag = 0.364 , 

04 = 0.299 , as = 0.281 , = 0.248 . 

This is again the case of a series with constant-sign coefficients, which is a difficult problem 
for resummation. For the approximant /I (2^), we find 

A21 = -0.990 , A22 = 0.778 , 7121 = -0.514 , 7122 = 0.164 . 

This gives the critical point Xc — |^2i|~^ = 1.010 and the critical index \n2i\ in a very good 
agreement with the exact values. In the next order, with an accuracy up to three decimal 
digits, we obtain 

^1 = -1 , As2 = 0.912 , ^3 = 0.363 , 

= -0.501 , 7132 = 0.091 , 7133 = 0.146 . 

The critical point Xc — \A3i\~^ — 1 coincides with the exact value, and the critical index |7i3i| 
is also practically equal to the exact index 1/2. The best Pade approximant, P[4/i](a:), is 
much less accurate, yielding the critical point Xc = 1.064 and a bad estimate for the critical 
index, equal to 1 instead of 1/2. The accuracy of the approximants f^ix) and P[4/i](a;) are 
compared in Fig. 10. 



4.3 Critical behavior in nonlinear differential equation 

Critical behavior may arise in solutions of differential equations. For instance, the Ruina- 
Dieterich law of solid friction, given by Eq. (27), is qualitatively different for 771 < 1 and 

m > 1. For m < 1, the solution has no zeros. But for m > 1, the solution becomes zero 
at a critical time approached with the critical index 1/m. Let m = 1.5 and all other 
parameters be the same as for the noncritical case considered in the previous section. Then 
the critical time is tc — 0.33004. For the short-time expansion in powers of time t, we have 
fo = /(O) = 0.5 and 

ai = —1.776 , a2 = —1.256 , = —1.706 , 

04 = -3.024 , as = -6.114 , as = -13.388 . 
The asymptotic series strongly diverges. Constructing the factor approximant /^C^), we get 

A21 = -3.049 , .422 = -0.558 , 7121 = 0.616 , 7122 = -0.181 . 
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This gives the critical time tc = \A2i\~^ = 0.328 and the critical index n2i, estimating rather 
well the corresponding exact values tc — 0.330 and 1/m = 0.667, respectively. For the 
parameters of /g (t) , we derive 

^31 = -3.036 , As2 = -1.676 , A33 = -0.225 , 

nsi = 0.630 , = -0.046 , n^s = -0.269 . 

The critical time tc = \Asi\~^ = 0.3294 approximates the exact numerical value with a 
good accuracy, the error being only —0.19%. The related critical index n^i is also close to 
the exact index 0.667. The most accurate Fade approximant P[3/3](t) yields a much worse 
approximation, with the critical time tc — 0.34134 and the critical index 1 instead of 2/3. 
The relative accuracies of the approximants f^ix) and P[2/2]{x) are demonstrated in Fig. 11. 



4.4 Two examples of critical phenomena in statistical physics 

An analysis of critical behavior would not be complete without considering critical phe- 
nomena of statistical mechanics. Consider, for instance, the three-dimensional spin-1/2 
Ising model with a simple cubic lattice. Thermodynamic characteristics of this model 
can be presented in the form of high-temperature expansions in powers of the parameter 
V = tanh( J/ZcbT), where J is the exchange integral and T, the temperature. For example, 
the second derivative of the susceptibility with respect to an external field 



^(2) _ d\ 

AO — 



(30) 



H=0 



is known [44] as the scries expansion containing the powers of v up to v^'^. The first coefficients 
of this expansion are /o = — 2 and 

ai = 24 , 02 = 318 , 03 = 3240 , 04 = 28158 , 

05 = 220680 , 0,6 = 1604406 , ar = 11029560 , ag = 72559422 . 

The coefficients are of constant sign and are growing fast. Applying the method of the factor 
approximants, we have in the second order 

A21 = -4.602 , A22 = 7.373 , risi = -4.300 , 7122 = 0.571 , 

where we limit ourselves by three decimal digits. This gives the critical point Vc = \A2i\~^ = 
0.217 at which the function diverges with the critical index |n2i|. These values can be com- 
pared with accurate Monte Carlo calculations [45] yielding the critical point Vc — 0.218092 
and with a refined analysis by means of integral approximants [44] giving the critical index 
4.37. In the third order, again with an accuracy up to three decimals, we obtain 

>l3i = -4.601 , A32 = 7.375 , ^3 = -44.934 , 

7131 = -4.301 , n32 = 0.571 , 7133 = . 
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Then the critical point is Vc — |A3i|^^ = 0.217 and the critical index is |n3i|, which are 
practically the same as in the second approximation. The best Fade approximant P[i/4](f) 
results in rather inaccurate values for the critical point = 0.153 and the critical index 
equal to 1. The fourth-order factor approximant does not lead to a noticeable change of 
the critical parameters, as compared to the third-order one, also giving Vc — 0.217 and the 
critical index 4.301. The percentage errors of those estimates are about 1%. 

As another example of high-temperature series expansions, let us consider such series for 
the (2 -|- l)-dimensional Ising model, defined by the Hamiltonian 

H-E{l-c7f)-xY^ala]-h^ a} , (31) 

i <ij> i 

in which erf, with a = 1, 2, 3, are Pauli matrices; the index i enumerates the sites on the two- 
dimensional spatial lattice; < ij > denotes nearest-neighbour pairs of sites; x is an effective 
coupling parameter corresponding to the inverse temperature in the Euclidean formulation 
of field theory; and h stands for magnetic field. Here, we consider the triangular lattice. Let 
us take, e.g., a series in powers of x for the mass gap F at zero magnetic field, 

F = Ei-Eo^ F{x) , (32) 

which is the difference between the energy of the first excited level and the ground-state 
energy. The series (2) for this case [46] has the coefficients fo — 2 and 

ai = —3 , a2 = —3 , — —5.25 , 04 = —15.75 , 

as = -49.265625 , = -173.3554688 , 07 = -627.602783 , og = -2397.718506 . 

For the parameters of the second-order factor approximant F2{x), we find 

A21 = -4.7404 , A22 = 1.7404 , riai = 0.6582 , 7122 = 0.0691 . 

The critical point, following from here, is Xc = |^2i|~^ = 0.2110. The mass gap tends to zero 
at as F ~ {xc — xY, with the critical index u = n2i. For the third-order approximant 
F^{x), we have 

A31 = -4.7826 , A32 = -3.5899 , A33 = 1.5136 , 

7131 = 0.6211 , 7132 = 0.0458 , 7133 = 0.0890 . 

Hence, the critical point is Xc = |^3i|~^ = 0.2091 and the critical index is u — 7131. In the 
next order, for the parameters of F^{x), we obtain 

^41^-4.7629, ^42 = -2.1880, ^43 = 3.9403, ^44 = 1-2670, 

7141 = 0.6445 , 7142 = 0.0375 , 7143 = 0.0007 , 7144 = 0.1178 . 

This results in the critical point Xc — |^4i|~^ = 0.2100 and the critical index u — 7141. These 
values can be compared with those summarized in Ref. [44], where Xc varies between 0.2097 
and 0.2098, while the index 1/ is between 0.627 and 0.641, which is very close to our results. 
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5 Conclusion 



By employing the self-similar approximation theory, we have derived a new class of approxi- 
mants, which, because of their form, we call the self-similar factor approximants. All control 
parameters of this class of approximants are defined by means of the accuracy-through-order 
relationship. These approximants reproduce exactly a wide class of functions from the sole 
knowledge of their asymptotic expansions. For other functions not from the exactly repro- 
ducible class, the factor approximants provide a very high accuracy, which is essentially 
higher than that given by the best Fade approximants constructed with the same number 
of asymptotic terms. 

The self-similar factor approximants are able to reproduce with a good accuracy various 
kinds of functions, diminishing and increasing, monotone and non-monotone, on finite or 
infinite intervals. The asymptotic series, used for constructing the factor approximants, 
can have their coefficients with alternating signs or with constant signs. The series can 
also be strongly divergent. Although the factor approximants are derived from asymptotic 
series for a variable in the vicinity of zero, they extrapolate well to the behavior of the 
sought functions at infinity. The approximants are able to predict the occurrence of critical 
phenomena, providing accurate values for both the critical points and critical indices. 

Note that a natural generalization involves combining the factor approximants with the 
self-similar root approximants and exponential approximants. The latter can be treated 
as a limiting case of the factor approximants, since the renormalized factor function (8) 
tends, as rik^ — > cxd, to an exponential. The self-similar exponential approximants, as has 
been recently shown [47], enjoy the property of exactly reconstructing exponential functions. 
Therefore, the factor approximants, together with their limiting exponential forms, and being 
combined with the self-similar root approximants, provide a powerful tool for an accurate 
reconstruction of a very wide class of functions. 

At first glance, it may appear mysterious that, knowing solely the behavior of a function 
in the vicinity of zero, it is possible to correctly predict its behavior at infinity or near a 
critical point. However, there is no miracle here. The coefficients of an asymptotic series 
contain a great deal of information about their parent function, provided that the latter are 
differentiable up to sufficiently high order. Then, the main problem is that this information is 
hidden, encoded. And one needs to possess a guide and a key for decoding that information. 
The idea of group self- similarity for subsequent approximations [17-27] serves as a guide 
pointing at the general properties of the sought function, which allows for the extrapolation 
of the given approximations. And the self-similar factor approximants is a practical key for 
realizing this extrapolation. 

Since the self-similar approximation method presented here is capable of capturing the 
features of rather complex functions with very good accuracy from the knowledge of a few 
numerical coefficients, one can view this approach as a complexity reduction scheme, or better 
as an encoding-decoding scheme. In the language of algorithmic complexity theory [Q, the 
differentiable functions with only isolated critical points have a low degree of complexity. 
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Figure Captions 



Fig. 1. Percentage errors of the self-similar factor approximant (solid line) and of 

the Pade approximant i-*[4/2](a;) (dotted line), as compared with function (13). 

Fig. 2. Percentage errors of the self-similar factor approximant f^{x) (solid line) and 
of the Pade approximants P[i/5](x) (dashed line) and Pi2/4]{x) (dotted line), compared with 
function (14). 

Fig. 3. Percentage errors of the factor approximant /|(a;) (solid line) and of the Pade 
approximants P[i/5]{x) (dotted line) and P[2/4]{x) (dashed line), compared with function (15). 

Fig. 4. Percentage errors of /|(a;) (solid line) and P[i/5](a;) (dashed line), approximating 
function (16). 

Fig. 5. Percentage errors of /3(a;) (solid line) and P[2/2]{x) (dotted line), compared with 
function (17). 

Fig. 6. Percentage errors of (solid line) and P[2/3]{9) (dashed line), compared with 
the numerical values of integral (18). 

Fig. 7. Percentage errors of the approximants for the expansion factor 0:3(2;) (solid line) 
and P[3/3]{z) (dashed line), compared with numerical values. 

Fig. 8. Percentage errors of the approximants for friction, f^{t) (solid hne), P[2/2]{t) 
(dotted line), and P[3/3](i) (dashed line), compared with an exact numerical solution of Eq. 
(27). 

Fig. 9. Percentage errors of fs{x) (solid hne) and P[3/2]{x) (dotted hne) as compared to 
function (28). 

Fig. 10. Percentage errors of f^{x) (solid hne) and P[4/ij{x) (dotted hne) with respect 
to function (29). 
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Fig. 11. Percentage errors of /^{x) (solid line) and P[s/3]{x) (dotted line), compared 
with the numerical solution of Eq. (27). 
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